Complete time-dependent treatment of a three-level system 
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Both unitary evolution and the effects of dissipation and decoherence for a general three-level 
system are of wide interest in quantum optics, molecular physics, and elsewhere. A previous paper 
presented a technique for solving the time-dependent operator equations involved but under certain 
restrictive conditions. We now extend our results to a general three-level system with arbitrary 
time- dependent Hamiltonians and Lindblad operators. Analytical handling of the SU(3) algebra 
of the eight operators involved leaves behind a set of coupled first-order differential equations for 
classical functions. Solution of this set gives a complete solution of the quantum problem, without 
having to invoke rotating-wave or other approximations. Numerical illustrations for multiphoton 
couplings and quantum control are given. 
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I. INTRODUCTION 

Three-level systems arise in many physical contexts. A 
spin-1 particle in an external magnetic field, three states 
of atoms or molecules coupled by lasers or other interac- 
tions, and oscillations among the three neutrino flavors 
are just a few varied examples that lead to a Schrodinger 
equation for the wave function or, equivalently, for the 
evolution operator U (t) , 

idU(t)/dt = H(t)U{t), U(0) =T. (1) 

Here, and throughout, we set K = 1. In a matrix repre- 
sentation, H(t) is a 3x3 matrix. As examples, population 
trapping and dispersion was considered through 

/ Gl \ 
H = Gi Ai G* 2 , (2) 
\ G 2 Ar + A 2 J 

where the coefficients G represent couplings between level 
one and two and between two and three while the diago- 
nal elements give the energy levels and detunings. Mul- 
tiphoton coupling in molecular systems has used Q 

/ -Eh. n 12 (t) o \ 
h= n 12 (t) o n 23 (t) , (3) 

V o n 23 (t) s J 

while recent treatments of quantum control in the 
rotating-wave approximation considered 0, Q 

/ o n la (t) o \ 
h = n 12 {t) a n 23 (t) . (4) 

V o n 23 {t) o J 

With constant coefficients for the diagonal entries, the 
two Hamiltonians in Eq. © and Eq. are essentially 
equivalent upon shifting the zero of the energy scale and 
suitable identification of the constants. 



Some of these studies introduced decay phenomeno- 
logically by considering complex values for the diagonal 
entries. A fuller treatment of dissipation and decoher- 
ence proceeds, however, by replacing the unitary evo- 
lution equation in Eq. by master equations for the 
density matrix p. A widely used class of such equations 
is the Liouville-von Neumann-Lindblad equation 0, E| , 

ip = [h, P ] + ([L k p, 4] + \ L k, p4i) 

k 

= [H, P] - { L t L ^P + P L l Lk - 2L xP L l) .(5) 

k 

where an over-dot denotes differentiation with respect to 
time, and the are operators in the system through 
which dissipation and decoherence are introduced. Even 
though this can result in non-unitary evolution, the form 
of the equation preserves Tr(p) and positivity of proba- 
bilities. A more mathematical discussion of such "super- 
operators" and "dynamical semigroups" is given in 0. 
Note that preservation of the trace so that there is no 
overall dissipation hinges on the coefficient 2 in the last 
term above, any other value leading to a change in the 
value of the trace with time. 

In recent papers 0, > we have developed a technique 
for solving Eq. JQ) and Eq. (JSJ). This method, which 
seems to have been independently rediscovered several 
times over the decades, with the earliest reference go- 
ing at least back to Wei and Norman separates the 
operator aspect from the time dependence by writing 

U(t)='[[exp[-i f i j (t)A j }, (6) 

j 

where the Aj are operators chosen to be time indepen- 
dent while the pj (t) are classical functions that carry all 
the time dependence. Our constructive procedure con- 
sists of inserting Eq. JBJ into Eq. whereupon the 
derivative on the left-hand side can be rearranged to 
have the same structure as the right-hand side through 



2 



repeated application of the Baker-Campbell-Hausdorff 
(BCH) identity and choosing the p,j(t) suitably. 
These functions are seen to satisfy a set of well-defined 
first-order differential equations. All the operator aspects 
are handled analytically in the BCH manipulations while 
the time dependence aspect is isolated into the equations 
for the n(t). These may require numerical solution but 
that is readily carried out through a simple MATHE- 
MATICA program |Dj . The same procedure also applies 
to the Lindblad structure of Eq. JSJ after first recasting 
it into the same form as Eq. namely, 



if,{t) = C(t) v (t), 



(7) 



where r)(t) are suitable linear combinations, (n 2 — 1) in 
number, of the density matrix elements themselves. Such 
a translation of Eq. JSJ into Eq. Q passes, of course, 
from the n x n matrix structure of H and Lk into (n 2 — 
1) x (n 2 — 1) matrices for the operators in C. Thus, the 
Lindblad equation for a two-level system is also of the 
3x3 form of Eq. (JTJ and Eq. Q Q. 



II. GENERAL SOLUTION OF A THREE-LEVEL 
PROBLEM 

The above construction requires that all the operators 
that arise as a result of the BCH application be included 
in the set of operators Aj in Eq. ©. Clearly, for the 
most general three-level problem, this requires eight lin- 
early independent matrices which, along with the unit 
matrix, provide a complete description. The eight trace- 
less matrices used in a standard description of SU(3) E2 
are a convenient choice. In our previous work we 
restricted our attention to forms of H and Lk for which 
a smaller subset, namely three, sufficed: 



-4, 




-4„ 




.4, 




(8) 



This set of three, displaying also their A notation in the 
SU(3) literature [13, coincide with the SO (3) angular 
momentum operators of three-dimensional rotations and 
close under mutual commutation between them. Thus, 
only three such terms and the unit operator are neces- 
sary in Eq. ©, affording a considerable simplification. 
(It helps to replace A x and A y by their linear combi- 
nations A x ± iA y .) In particular, the resulting set of 
equations for the three /j, are readily amenable to solu- 
tion , consisting of a Riccati equation (first-order and 
quadratically nonlinear) for one of them whose solution 
then leads to simple quadrature solutions for the other 
two. The H and Lk considered in [y, |8j permitted such 



a use of the SO(3) sub-group of SU(3) to simplify the 
solution of Eqs. QJ and JSJ. 

However, Hamiltonians such as those in Eqs. @-(QJ 
cannot be expressed in terms of just the three operators 
in Eq. (JSJ. Similarly, a two- level Lindblad equation ^1 
with H = ^e(t)a z + J(t)a x , L = VTct z , gives an Eq. (J7J) 
of the form |(j 



dt 
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(9) 



with Ei = E2 = —iT, Again, the C operator requires use 
of the full set of eight operators of the SU(3) algebra. 

Such a set of eight includes, besides the three in Eq. JSJl, 
three other off-diagonal matrices [l2|. 
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(10) 



and two diagonal ones, 




1 



1 



— 010 

v3 \ -2 



(11) 



All eight operators and the unit operator are needed 
in Eq. © for the general solution. Once again, the lin- 
ear combinations that simplify the commutators and, 
therefore, the final set of equations for are a± = 
±i\ 7 ),b± = i(A x T «A 2 ),c± = i(A 4 ±iA 5 ), and 
0-3 = 5(v3As — A3),c 3 = |(v3Ag + A3). Table 1 shows 
the commutators between this set of eight operators. The 
three sets of triplets {a, b, c} coincide to within factors 
ofj with the sets {V, T, U} used in the SU(3) literature 

EI 

Flexibility in the order of the various factors in Eq. © 
is a feature of our technique. Inspection of Table 1 per- 
mits optimal ordering such that each application of the 
BCH identity generates at most two terms. We choose, 
thereby, the order 



U(t) 



& -iS e -in s b+ e -iH7b- e ~ifj, 6 c + e ~ip, 5 c^ 

x g— ipsa+g— ip,2d- e -ima 3 e -ip.AC 3 ^2) 



Evaluating idU/dt, re-arranging through use of the 
BCH identity into the form of an operator sum 
acting from the left on U(t), the coefficients of 
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TABLE I: Table of commutators. With operators Oi in the first column and Oj in the top row, each entry provides the 
commutator [Oi,Oj]. 



b + , c+, c_, o+, a_, 03, C3 in the operator sum are, re- 
spectively, 

/i 8 + »M5S + [4,w + z/i 8 (-i/i 3 r + /ii -2/i 5 /i 7 s -/i 4 + «M6«) 
w = fir — ifJ-er + i/is/iyS —ifi7(—ifisr + (i\ + i/i 6 ?; — /i 4 ) 
u = M6 + /^^ + ipe{-i^r + fii + 2/i 4 ) -1/17(1 -^^s 

v + i/i, 8 (r + i/^-u) 

-i/i 8 u + (1 - /i5M6)s 

r + i/i7U 

/ii - z/i 8 w - «/i 3 r - /i 5 M7S 

/i 4 + z/i 8 ui - l^JLQV + /i 5 /i 7 S, (13) 

along with J multiplying X. We have denned for conve- 
nience, 

r = (12 - «/i2(M4 + 2/ii), 

s = /i 3 + //3V + i/i 3 (/i 4 + 2/ii) 

v = /i 5 - i^ 5 (-ifi 3 r + /i 4 + 2/i 4 ). (14) 

The above set of expressions can be matched to any H 
in Eq. or £ in Eq. J7J with nine arbitrary time- 
dependent entries to provide defining equations for the 
[i and (5 in Eq. (|l"2"|> . This completes the general solution. 

Thus, for the matrix in Eq. @, which is a linear com- 
bination of (a±, 6±, 0,3, c 3 ), we have 

/i 8 + Ms 6 + i^s(E 2 -E 1 ) =-e- i/z 5 m, 
/i 7 -ifi 5 f4m -2/i 7 /i 8 e - £i 7 (-E2 - -Ei) -«>6" =-e, 
/i 6 + i/x 6 /i 8 (2/x 6 J + ie)-i^ 5 ^ 6 fj, 7 m+ i^Ei = 2i/i 7 J, 
/i 5 + i^l^jm -?^ 5 /i 8 (4/i 6 J + ie) -i^Ei = -2ifi s J, 

/i 4 =2/i 6 /i 8 J + 2/i 8 e-/i 5 /i 7 m + ^(2^1 --E2), 

A3-A*3 n + ^3(^*5^7™ + M8[2/i6"/-*e] + ^2) = 
/i 2 -i^2M5M7W -i/i 2 /i 8 (2/i 6 J ~«e) + 2/x 2 /x 3 n -i(i 2 E 2 = n, 

(ix= i/i 3 n +M5j"7W-i^ 8 e +i(2£ 2 -£a), (15) 

with 5 = (E^ + £ , 2 )/3, where we have defined m = 
2J/(1 — /j,5/j,e),n = 2J(1 — /X7/i 8 ). The above eight equa- 
tions fall into two groups, the first four involving only 



fi$ — fig. They may be solved first and then serve as 
inputs for solving the remaining four equations. This 
structure was anticipated in our previous work ||, be- 
cause the 8x8 matrices in Eq. Q) that follow from the 
operators (a±, 6±, 03, C3) are in 4x 4 block-diagonal form. 
The above set of equations are, of course, more com- 
plicated than in [y, |8j, being more nonlinear but note 
again the characteristic structure of no higher power than 
quadratic of any \i. 

After completion of our work, we became aware of 
closely related papers that use the Wei-Norman Q pro- 
cedure for three-level problems 0, [^. In particular, 
|15| uses a product of nine exponentials involving the op- 
erators {T, U, V} of SU(3) with an arrangement slightly 
different in order from ours in Eq. i|12|) . Our choice seems 
to yield a less complicated set of equations in Eq. i|15|) 
that define the fi functions. 



III. LINDBLAD EQUATION FOR A 
TWO-LEVEL SYSTEM 

The most general two-level problem with an arbitrary 
time-dependent Hamiltonian and any choice for the Lk 
in the master equation in Eq. (jSJ leads to a 3 x 3 equa- 
tion for the density matrix elements such as Eq. 10 and 
is thereby solved through the set in Eq. i|15[l . Fig. 1 
provides an illustration for particular choices of the pa- 
rameters in Eq. ©. The left-hand side of the panels 
show perfect agreement with the solutions obtained by 
a different numerical procedure in [T^ . The right-hand 
side shows for comparison our earlier results [g where 
the decoherence part involving T was simplified so as to 
reduce the second term on the right-hand side of Eq. © 
to a unit-diagonal contribution on the right-hand sides of 
Eq. (0) and Eq. (0. As expected, the difference between 
the two lies mainly in the off-diagonal density matrix el- 
ement and the rate of damping. The entropy, while still 
rising monotonically to In 2 as the system evolves from a 
pure state to a mixed one, also differs in these two mod- 
els of the decoherence. The comparison suggests that the 
simplified model for decoherence introduced in 6] works 
reasonably well. Since it is much easier to implement, 
especially for higher n, requiring only n-dimensional ma- 
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FIG. 1: Diagonal and off-diagonal elements of the density 
matrix and entropy S for a two-level system described by 
Eq. © and Eq. Q. With e(t) = Acos(ujt),2J = Scos(Ht + 
5), results are shown in the left-hand panels for an initial 
pure state 1, and A = 45, B = 6,w = LO = 0,5 = 0,T = 
0.3. These results coincide with those of The right-hand 
panels show for comparison the results of a simplified model 
for decoherence in 



trices rather than n 2 — 1, we will so use it below for 
three-level systems. 



IV. APPLICATIONS TO THREE-LEVEL 
SYSTEMS 



In a study of population trapping and dispersion in 
three-level systems, pj considered the Hamiltonian in 
Eq. (J2J for both "lambda" and "ladder" systems with 
Ai = L02 — u>i — u) a , A2 = 0J3 — L02 i u>b, the signs cor- 
responding to lambda and ladder arrangements, respec- 
tively, of the three levels 1,2, and 3. Solutions were 
developed for special cases such as equal amplitudes of 
the driving fields and a phenomenologically introduced 
damping through a negative imaginary piece in u>2- Fig. 
2 presents our results through Eq. JSJ and Eq. (|15|) with 
the decay introduced as per the model discussed above. 
In this model, Eq. JSJ) reduces to Eq. @ with the de- 
cay occurring as — ir times the unit operator in C The 
results coincide with those of .1]. But our procedure ex- 
tends readily to arbitrary values of the parameters in the 
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FIG. 2: Populations of a three-level system described by 
Eq. starting with initial state 1. The parameters are 
|Gi| = |G 2 | =0.5 and (a) A x = 0.5 - O.Oli, A 2 = 0.5 + O.Oli, 
(b) Ai = —A2 = 5 — i. The intermediate level 2, shown by a 
thick line, remains essentially at zero for all time in (b). The 
results coincide with those of Q. 



Hamiltonian and can be applied also to varying ampli- 
tudes, frequencies, and phases of the two driving fields. 
The same numerical solutions of Eq. I|15[l , when inserted 
into Eq. (|12Jl provide a full solution of the time evolution 
of all density matrix elements, whatever the time depen- 
dences and values of the parameters in the Hamiltonian 
in Eq. (0) . Note in the bottom panel of Fig. 2 that the 
intermediate level 2 has negligible population through- 
out, the population flopping back and forth between the 
extreme levels. 

Another more recent paper [jj presented an analyti- 
cal model of three-level systems with the Hamiltonian in 
Eq. (0 but again under restrictive conditions where the 
two driving fields have the same hyperbolic-secant time 
dependence in at, differing only in amplitude, and with 
the real part of the detuning A set equal to zero. The 
value of a sets the time scale. In Fig. 3 we present our 
results for a wider range of parameters. Our results for 
the diagonal components (the off-diagonal ones are not 
displayed but are also available in our calculations) of the 
density matrix are presented as functions of time but are 
equivalent to the display in terms of detuning given in 
|3| . The change in form with increasing values of the pa- 
rameters a — yJ~A 



A\ z a, S — hAcr, and 7 



iTer are 

interesting. The population of states 1 and 3 at large t 
depend critically on these parameters. Here the A are the 
amplitudes of the two driving fields, and Y the damping 
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FIG. 3: Populations of a three-level system described by 
Eq. starting with initial state 1. Thick and thin lines 
describe levels 2 and 3, respectively, and a dashed line the 
population of level 1. Values of parameters described in the 
text are (a) a = 2, 7 = 0, S = 0, (b) a = 2, 7 = 0, 5 = 0.866, 
(c) a = 5, 7 = 0, 5 = 0, (d) a = 5, 7 = 0.5, 8 = 0. 



(with -«r added to A in Eq. |@J)) as in la. Our results 
provide a method for exploring a broad range of param- 
eter values and time dependences in Eq. fllf. 



Finally, as yet another application, we consider a 
very recent paper pj that used a Hamiltonian encod- 
ing scheme based on computational control techniques 
for handling the Dyson series in time-dependent evolu- 
tion. The results presented were for zero detuning A in 
which case the problem in Eq. 0} actually reduces to the 
much simpler 3x3 problem we considered previously • 
But, we present in Fig. 4 much more general results for 
a wide range of parameters in that Hamiltonian. Decay 
is again introduced through the simplified model (with 
— iT added to the diagonal entries in Eq. (J2J) that per- 
mits application of the 8x8 problem in Eq. and 
Eq. (JT3J. The form of the driving fields in Eq. J2J is 
^12 = (j4i /2) exp(— [(t — ti)/a] 2 ) and a similar expres- 
sion for Q23 with subscript 2 for the parameters. Starting 
with initial population in the state 1, Fig. 4 shows the 
subsequent evolution of the three populations. Larger 
values of a show initial oscillations as the population is 
transferred from 1 to 3 asymptotically. The effects of 
damping when all three populations equalize, as well as 
the effect of increasing A are also shown. 
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FIG. 4: Populations of a three-level system described by 
Eq. , starting with initial state 1. Values of parameters 
described in the text are A2 — A\ = 2.5, ti = 12, ti — tj — a, 
and (a) a = 8, (b) a = 3, (c) a = 1.5, (d) a = 0.9, all 
with 7 = and A = 0, (e) a = 3,7 = 0.15, A = 0, (f) 
a = 3,7 = 0, A = 1.5. 



SUMMARY 



We have developed a complete solution of the time evo- 
lution of three-level systems. With individual, arbitrary 
time-dependent entries in the 3x3 Hamiltonian, equat- 
ing those entries to the set in Eq. (|T3)l provides coupled, 
first-order differential equations for classical functions u. 
They are amenable to simple integration through pT| 
and the solutions, when inserted into Eq. I|12fl . provide a 
complete solution of the quantal problem. Thereby, there 
is no need to invoke any familiar approximation schemes 
used such as the rotating-wave approximation. Master 
equations for two-level systems with dissipation and de- 
coherence, when the most general form is of 3 x 3 matrix 
form, are also solved in the same manner. A simplified 
model of decoherence allows solution of three-level prob- 
lems as well. Illustrations are given of a variety of appli- 
cations from the recent literature for population trans- 
fer and trapping in quantum optics and for multiphoton 
transitions in molecules. 

This work has been supported by the U.S. Department 
of Energy under Grant No. DE-FG02-02ER46018. 
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